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Abstract 

We calculate the structural and electronic properties of OsB2 using density functional theory 
with or without taking into account spin-orbit (SO) interaction. Our results show that the bulk 
modulus with and without SO interaction are 364 and 365 Gpa respectively, both are in good 
agreement with experiment (365-395 Gpa). The evidence of covalent bonding of Os-B, which plays 
an important role to form a hard material, is indicated both in charge density, atoms in molecules 
analysis, and density of states analysis. The good metallicity and hardness of OsB2 might suggest 
its potential application as hard conductors. 
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Hard materials are of great interest due to their superior properties of higher compres- 
sional strength, thermal conductivity, refractive index, and chemical stability beside higher 
hardness.— Besides the two well-known super-hard materials of diamond and cubic boron 
nitride (c-BN), experimental and theoretical efforts are devoted to searching a new class 
of hard meterials. Covalent materials are much better candidates for high hardness than 
ionic compounds because electrostatic interactions are omnidirectional and yield low bond- 
bending force constants, which result in low shear modulus?^ Some potential hard materials 
such as C3N4 polymorphs,— Si3N4,— BC2N,— and B60,— have been proposed and investigated 
intensively in recent years. 

The pure Osmium (Os) has an exceptionally high bulk modulus (395-462 GPa) while 
its hardness is of only 400 kg/mm^.— The high bulk modulus of Os mainly due to the 
high valence electron density, while the low hardness is related to the metallic bonds and 
the hexagonal close-packed (HCP) crystal structure of Os. It has been found that many 
transition metals, such as W, Ru, and Zr, which are soft in their pure metals, can be 
converted into hard materials by combining with small covalent bond-forming atoms such 
as boron, carbon, oxygen, or nitrogen . For Os related materials, Os02 has been 
predicted to be very stiff,^^ and some compounds like OsBq.h, OsB2, and OS2B3 have already 
been synthesized.— Recently, Cumberland et al. studied experimentally the mechanical 
properties of the OSB2 compound, and showed that it is an ultra-incompressible (365-395 
GPa) and hard (>2000kg/mm^) material.-^ 

In this brief report, we carry out density functional calculations on the electronic proper- 
ties and the bulk modulus of OSB2, using the highly accurate full potential linearized aug- 
mented plane wave plus local orbital (APW + LO) method as implemented in the Wien2k 
code.^^ The nonoverlapping muffin-tin (MT) sphere radiis of 2.02 and 1.57 bohr are used 
for the Os and B atoms, respectively. Expansion in spherical harmonics for the radial wave 
functions are taken up to 1=10. The value of YiMT^max (the smallest muffin tin radius 
multiplied by the maximum k value in the expansion of plane waves in the basis set) which 
determines the accuracy of the basis set used is set to 7.0. The total Brillouin zones are sam- 
pled with 6000 k points, i.e., 672 k points in the irreducible wedge of the Brillouin zone. For 
the exchange-correlation energy functional, we use the local density approximation (LDA) 
in our calculations. Since Os is a heavy element, a fully relativistic calculation is performed 
for core states, whereas the valence states are treated in both scalar and fully relativistic 
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fashions. The spin-orbit couphng (SOC) is self-consistently added via the second variational 
step scheme with relativistic pi/2 local orbitals included.— We find there is no spontaneous 
spin polarization in OSB2 in the framework of the local spin density approximation (LSDA) 
or LSDA+SO methods. Hereafter for the sake of simplicity, we only discussed the results 
from the spin-restricted calculations. 

OsB2 has an orthorhombic lattice (space group Pmmn, No. 59) with the experimental 
lattice parameters a = 4.6832 ± 0.0001, b = 2.8717 ± 0.0002, and c = 4.0761 ± 0.0001 A. In 
the orthorhombic structure, two Os atoms occupy the 2a Wyckoff sites and four B atoms 
occupy the 4f sites, as shown in Fig.^ We optimize both lattice and internal coordinates to 
get a relaxed structure for OSB2. The total energy per unit cell as a function of volume is 
calculated and the results are fitted using the Murnaghan equation of state to get the bulk 
modulus.— The corresponding equilibrium lattice parameters, bulk modulus from both LDA 
and LDA+SO methods are compared with the experimental result, as shown in Table lU We 
can see that the equilibrium lattice parameters calculated both in the absence and in the 
presence of SO are very close with the experimental value. And the bulk modulus calculated 
in the absence and in the presence of SO are 365 and 364 Gpa, respectively, both are in 
good agreement with the experimental result (365-395 GPa).^^ 

As discussed above, the hardness of a material is largely determined by the bonding 
type in the system. Thus the electronic properties of OSB2 such as the electron density, 
band structure, and the density of state (DOS) are calculated in the absence and in the 
presence of SO at the equilibrium geometries. We find that the electron density from the 
LDA calculation is almost the same as that from the LDA+SO calculation. Here we show 
the electron density for OSB2 from the LDA+SO calculation in Fig. |21 Two charge density 
contour plots (one in the plane including Osl, Os2, B2, and B3, and the other in the 
plane including Osl, Bl, and B2) are shown in Fig. |21^a) and (b), respectively. There are 
some electrons between Os atoms and theirs neighbor B atoms, indicating covalent bonding 
in OsB2. To gain a more detailed insight into the bonding characters of OSB2, we plot the 
difference density (the difference between total density and superposition of atomic densities) 
in Fig. |21^c). We can see that some electrons are transferred from both B and Os atoms to 
B-B bonds and Os-B bonds. Clearly, two neighbor B atoms form very strong covalent bond. 
And the bonding between the Os atom and its neighbor B atoms is also covalent. The 
fact that Os-B bond is covalent rather than ionic is also evident from Bader's atoms in 
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molecules (AIM) theory,— which indicates there is only about 0.05 electron transferred from 
an Os atom to B atoms. Moreover, there is some region where the charge density is rather 
low, implying the bonding in this system is directional but not metallic, as can be seen 
from Fig. |2^a) and (b). Since there are eight B atoms around an Os atom, forming a three 
dimensional network instead of a plane, the anisotropy in OsB2 is not expected to be very 
large. As discussed above, the strong covalent bonding of Os-B plays an important role to 
form a hard material, since the highly directional bonding is needed to withstand both elastic 
and plastic deformations. OsB2 might be potential superhard material for effectively cutting 
ferrous metals, including steel, since diamond reacts with steel producing iron carbide.— A 
recent report showed that the hardness of the Fe-based alloys is less than 1300kg/mm^,— 
which is far less than that of the OsB2(>2000kg/mm^).— 

Although the LDA electron density is almost indiscernible from that from the LDA+SO 
calculation, we note that there are some differences in the band structures and DOS between 
two results. Fig. a) and (b) show the band structures near the Fermi energy in the absence 
or presence of SO, respectively. Due to the spin-orbit interaction, some bands are split: some 
bands crossing each other in Fig. Efa) no longer cross in Fig. Efb), and more noticeably 
almost all bands in the direction of M-N-K are split in the LDA+SO result. Except the 
band splitting in the presence of SO, the both band structures have main features in common. 
Both Fig. Efa) and (b) show that the bands are dispersive but not fiat indicating not large 
anisotropy in OsB2. We note that there are several bands crossing the Fermi level which are 
rather dispersive in the whole Brillouin zone, indicating good electronic mobility in OsB2. 
The metallicity is uncommon in superhard materials. Most of the superhard materials 
are insulators or semicondunctor.— This special character of OsB2 should bring it special 
application on electron conductivity. Fig. |3]show the total and partial DOS with or without 
SOC which are obtained using a modified tetrahedron method of Blochl et al.M' We can 
see some minor differences between the results with or without SOC. Fig. Wi^) also implies 
the good metallicity due to the substaintially large total DOS at the Fermi level. From 
the calculated partial DOS in Fig. |3J one can see that the electronic structure of OsB2 is 
governed by a strong hybridization between the Os-d and B-p states, while with a rather 
small contribution from the Os-p and B-s states. This strong hybridization of Os-d and B-p 
also indicates the strong covalent bonding of Os-B. 

Besides the bulk modulus, we want to mention that the shear modulus is also relevant 
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for the non-cubic OSB2 materials. The shear modulus can be computed from the elastic 
constants. Unfortunately, the force calculation when including the spin-orbit coupling has 
not been implemented in the Wien2k package yet. As an approximation, in the calculation 
of elastic constants, we keep the relative positions of the atoms fixed at small elastic strains 
of the lattice. Both calculations without SOC or with SOC are carried out. The procedure 
for computing elastic constants is similar to that used by Zhou et al.."^"^ For comparison, we 
also compute the elastic constants using VASP^ with LDA but no SOC included. In these 
calculations, we perform two series of calculations with atomic internal coordinates fixed or 
relaxed, respectively. Our results reported in Table HH We can see that the Voigt shear 
modulus G calculated using the Wien2k package (235.1 and 251.8 GPa when including SOC 
or not, respectively.) are in consistent with that (253.6 GPa) calculated using VASP with 
fixed atomic internal coordinates, but are much larger than that (197.1 GPa) calculated 
using VASP with a relaxation of the atomic internal coordinates. The results produced by 
VASP are in good agreement with those reported by Zhou et al.M Thus we find that is not a 
very good approximate that computing the shear modulus with atomic internal coordinates 
fixed. The accurate shear modulus with SOC included might be computed in our future 
works. 

In conclusion, the mechanical and electronic properties of OSB2 are studied using the 
highly accurate full potential linearized APW + LO method. The OSB2 is found to be 
metallic, which is uncommon in hard materials. The calculated bulk modulus with and 
without SOC are 364 and 365 GPa respectively, which are in good agreement with the 
experiment (365-395GPa).— The evidence of covalent bonding of Os-B, which plays an 
important role to form a hard material, is manifested both in the difference charge density, 
AIM analysis, and DOS analysis. Our results indicate that OSB2 might serve as hard 
conductors. 
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TABLE I: Structural parameters and bulk modulus of OsB2. 



LDA LDA+SO Experiment 



Lattice parameters (A): 

a 4.6433 4.6581 4.6832±0.000li4 

b 2.8467 2.8700 2.8717±0.0002ii 

c 4.0432 4.0560 4.0761±0.000lii 

Bulk modulus(GPa) 365 364 365-395^^ 



TABLE IL Calculated elastic constants and Voigt shear modulus G of orthorhombic OsB2. Unless 
otherwise stated, all results are computed with atomic internal coordinates unrelaxed. All values 
are in units of GPa. 



Wien2K LDA Wien2K LDA+SO VASP LDA VASP LDA (relaxed) Other VASP LDA (relaxed)22, 



Cii 


628.9 


608.5 


622.4 


597.2 


597.0 


Cl2 


194.7 


198.3 


185.8 


188.7 


198.1 


Cl3 


235.0 


220.9 


225.3 


217.5 


206.1 


C22 


627.8 


590.3 


624.9 


584.5 


581.2 


C23 


126.8 


129.9 


124.1 


164.0 


142.6 


C33 


923.2 


855.7 


910.4 


833.8 


825.0 


C44 


185.5 


175.1 


182.0 


80.2 


70.1 


C55 


313.5 


292.6 


311.4 


214.5 


212.0 


Cm 


218.6 


205.9 


233.9 


209.2 


201.3 


G 


251.8 


235.1 


253.6 


197.1 


193.8 
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FIG. 1: (Color online) The crystal structure for OsB2. The lattice vectors are denoted as a, b, 
and c. 




FIG. 2: (Color online) (a) The valence charge density contour (in e/a.u.^) plot in the plane including 
Osl, Os2, B2, and B3. (b) and (c) show respectively the valence charge density and difference charge 
density contour (in e/a.u.^) plot in the plane including Osl, Bl, and B2. The charge density shown 
here is the result from the LDA+SO calculation. The LDA charge density has similar character. 
Osl, Os2, Bl, B2, and B3 refer to the atoms labeled in Fig. ^a). 
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FIG. 3: (a) Band structure of OsB2 from the LDA calculation, (b) Band structure of OsB2 
from the LDA+SO calculation. T = (0,0,0), M = (0,0.5,0), = (0.5,0.5,0), K = (0.5,0,0), 
A = (0,0,0.5). 
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FIG. 4: (Color online) Partial and total DOS of OsB2 from the LDA and LDA+SO calculations. 




